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Abstract. We present a sample of 8 nearby Seyfert 2 galaxies observed by HST and Chandra. All of the sources present soft 
X-ray emission which is coincident in extension and overall morphology with the [O III] emission. The spectral analysis reveals 
that the soft X-ray emission of all the objects is likely to be dominated by a photoionized gas. This is strongly supported by 
the 190 ks combined XMM-A'evvton/RGS spectrum of Mrk 3, which different diagnostic tools confirm as being produced in a 
gas in photoionization equilibrium with an important contribution from resonant scattering. We tested with the code CLOUDY 
a simple scenario where the same gas photoionized by the nuclear continuum produces both the soft X-ray and the [O III] 
emission. Solutions satisfying the observed ratio between the two components exist, and require the density to decrease with 
radius roughly like r^^, similarly to what often found for the Narrow Line Region. 

Key words, galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

Highly obscured objects are ideal laboratories to study the 
properties of the gas in the environment of Active Galactic 
Nuclei (AGN). The high energy X-ray spectrum is dominated 
by a heavily absorbed power law and, especially when the col- 
umn density is very large (as for the Compton-thick objects, 
withA^H > o"/ ' - 1 -5 xlO^'* cm"-), by Compton reflection and a 
strong iron fluorescence line. Both the last two components are 
likely to be produced in a compact (parsec-scale), Compton- 
thick and low ionized material, traditionally identified with the 
absorber along the line of sight (the so-called 'torus' envisaged 
in the Seyfert unification scenarios: lAntonuccil 1 993l) . 

On the other hand, the soft X-ray spectrum had to await 
the high resolution spectrometers aboard XMM-Newton and 
Chandra to reveal its real nature, at least in the three bright- 
est sources, where such an experiment is possible: NGC 1068 
(iKinkhabwala et al.1 l2002t iBrinkman et al.. .2002: Ogle et al 



2003 )■ Circinus dSarnbruna et alJ2001bl) and Mrk 3 dSako et al 



200ol iBianchi et al.1 l2005bt IPounds & Pa ge" 200?). In these 
cases, the spectrum is dominated by emission lines, mainly 
from He- and H-like K transitions of light metals and L tran- 
sitions of Fe. The 'continuum' observed in lower resolution 
CCD spectra was mainly due to the blending of these features. 
Different diagnostic tests agree that in all these sources the ob- 
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served lines are more easily explained if produced in a gas pho- 
toionized by the AGN, rather than in a hot gas in collisional 
equilibrium. 

The unprecedented high spatial resolution of Chandra 
added another key ingredient to the puzzle. The same three 
sources revealed that their soft X-ray emission was actually 
extended on hundreds of p c, thus being clear l y produced in a 
gas well beyond th e torus ('Y oung et alJl200ll ISambruna et alJ 
12001 at ISakoet alj EoOO). Moreover, the dimension and mor- 
phology of this emission closely resemble that of the Narrow 
Line Region (NLR), as mapped by the [O III] /15007 emission 
line. Since the NLR is also generally believed to be mainly con- 
stituted of gas photoionized by the nuclear continuum, this cor- 
relation deserves a deeper investigation. Some positive results 
in this d irection have already been published on sin gle sources 
(see e.g. lKinkhabwalaet alJ200llOgle et alJ2003l) . 

We present in this paper a sample of 8 Seyfert 2 galaxies 
with extended [Olll] emission in their HST images, which 
were also observed by Chandra (Sect.|2li. We show in Sect. 13. II 
that all of them present a soft X-ray emission with extension 
and morphology highly correlated with that of the NLR. 
Moreover, the spectral analysis generally suggests that the 
most likely origin for this emission is from a photoionized gas 
(Sect. 13. 2> . Therefore, we produce photoionization models to 
find out if the same gas can be actually responsible for both 
components (Sect. 13, and finally discuss on the results (Sect. 
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The cosmological parameters used throughout this paper 
are Hq = 70 km s"' Mpc"', Ap = 73 aii d = (i.e. the 
defauhones in XSPEC 12.2.0: lArnaudll996l) . 

2. Observations 

2.1. The sample 

Table[npresents the sources included in the sample under anal- 
ysis, along with some basic properties and the details on the 
Chandra and HST observations analyzed in this paper. The 
sample consists of al l the Seyfert 2 galaxies included in the 
ISchmitt et alJ ('2003a') catalog, with a Chandra observation. 
The complete Schmitt et al. (2003a) catalog is selected on the 
basis of 60 /im luminosities and warm 25/60 /im colors, in order 
to avoid starburst galaxies and orientation biases, and includes 
only Seyfert galaxies with z < 0.031. We excluded NGC 1068, 
because its Chandra image and co rrelation with optica l and ra- 
dio data were treated in detail by 'Youna et al.' ("2001*). On the 
other hand, we added another source, NGC 5643, whose [O III] 
image was discussed bv lSimpson etalJ(ll997h . 

2.2. Data reduction 
2.2.1. Chandra 

We report here for the first time on the Chandra ACIS obser- 
vations of all the sources i n the sample , with the exceptions 
of Mrk 3 dSako et alJl2000|). NGC 4388 (llwasawa et aDl2003l) 
and NGC 4507 (Matt et al. 2004). All data were reduced with 
the Chandra Interactive Analysis of Observations (CIAO) 3.2.1 
and the Chandra Calibration Database (CALDB) 3.0.3, adopt- 
ing standard procedures. In particular, a new evt2 file was cre- 
ated with acis_process_event s, adopting an observation- 
specific bad pixel file and dedicated background-cleaning pro- 
cedures in the cases of VFAINT observations. Moreover, filter- 
ing of periods of anomalous background levels and correction 
for ACIS contamination at short wavelengths were performed. 
Spectra were extracted with PSEXTRACT, adopting three types 
of regions for each source: 1) circular regions with radius of 
1 arcsec (the 'nucleus' hereafter); 2) annuli with inner radius 
of 1 arcsec and a different outer radius for each source, de- 
pending on the size of the extended region ('extended region' 
hereafter); 3) circular regions encompassing both. All images 
were corrected for any known offset, to achieve the best possi- 
ble accuracy in astrometry. 

All spectra were binned in order to have no less than 25 
counts in each spectral channel, allowing us to use the;\f^ statis- 
tics. On the other hand, 'local fits' were also performed in the 
5.2 5-7.25 keV e nergy range with the unbinned spectra, using 
the lCashI statistics, in order to better assess the pres- 

ence of the iron lines and measure their properties. We refer 
the reader to Gu ainazzi et alJ ( l2005bl) for details on this kind of 
analysis. 

In the following, we will always call 'soft' and 'hard' the 
X-ray spectra below and above 3 keV, respectively. Fit inter- 



vals depend on the quality of the spectra, but are in any case re- 
stricted to the range 0.3-9.5 keV, which is approximately where 
the ACIS CCDs are well calibrated. 

2.2.2. HST 

We downloaded the narrowband optical images from the 
MAST (multi-mission archive at STScI). The observations 
were made as part of different programs, and with different 
instruments, filters and exposure times. The log of HST ob- 
servations is reported in Tabled The images were processed 
through the standard OTFR (on-the-fly reprocessing) calibra- 
tion pipeline which performs analog-to-digital conversion, bad 
pixel masking, bias and dark subtraction, flat field correction 
and photometric calibration. The data, usually a set of two im- 
ages for each object, are then combined to reject cosmic rays. 
The images taken with the FOC are corrected for geometric 
distorsion by the pipelin e, while for WF PC2 data we used the 
multidrizzle script JKoekemoer et al..,2002) . 

A crucial issue is clearly represented by the relative as- 
trometry of the images from Chandra and HST. Chandra has 
a nominal position accuracy of 0.6 arcsec (at the 90% confi- 
dence level), while the absolute astrometry of HST is accurate 
to 1-2 arcsec. Since the NLR are extended on scales of few arc- 
sec, the systematic offset between the images can be relevant. 
Unfortunately, in none of the datasets included in our sam- 
ple point-like sources are present in the fields of both instru- 
ments that could help to align the two astrometries. We there- 
fore decided to match the two pixels with the largest count- 
rate, which should represent the nucleus in both wavebands. 
The uncertainty relating to the fact that these sources are ob- 
scured, which, in principle, would displace the optical from 
the X-ray peak, should not be very relevant, because the phys- 
ical size of the Chandra pixel (around 0.5 arcsec) is much 
larger than the scale where the obscuration takes place (see 
Table 0. The reliability of this method, at least for the aims 
of this paper, is supported by the case of NGC 1068, whose 
HST and Chandra peaks, independen tly align ed, are coinci- 
dent within the Chandra errors jYounget alJi200ll) . Mrk 3 
represents the only exception, since, as already noted by e.g. 
Ruiz et al. (2005), it lacks a bright central [Olll] knot. In this 
case, we assumed that the optical nucleus were in the center of 
symmetry of the observed structure, which also coincides with 
the pixel with the lowest count-rate. 

2.2.3. Xmm-Newton 

In this paper we present the high-resolution combined spec- 
trum of Mrk 3, because this is the only source in the sam- 
ple whose RGS spectrum has a quality good enough to pro- 
vide us with useful gas diagnostics constraints. Mrk 3 was ob- 
served 10 times by XMM-Newton between 2000 and 2002. 
In all these observations, the Radiation Grating Spectrometer 
(RGS ; Iden Herder et"ai]l200ll) was operating in standard spec- 
troscopic mode . Part of the data o f these obs ervations were dis - 
cussed by Bianchi et alJ (l2005bl) and .Pounds & Paed (l2005l) . 
We have reduced the data of all flie RGS XMM-Newton ob- 
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Table 1. The log of all the observations analysed in this paper. 



Name 


z 


1"" 


Chandra 






HST 




Fo.5-2keV 


F[oiii] 


[Oiii]/ 






(pc) 


Date 


Instr.* 


Exp.' 


Prog. ID 


w 


Exp."^ 


(10-" cgs) 


(10-" cgs) 


softX 


NGC 1386 


0.0029 


60 


2003-11-19 


A-S 


20 


6419 


W2/F5 


800 


1.8 


5.1 


2.8 


Mrk3 


0.0135 


275 


2000-03-18 


A-S/H 


100 


5140 


F/F5 


750 


4.7 


20.5 


4.4 


NGC 3393 


0.0125 


255 


2004-02-28 


A-S 


29 


3981 


W/F5 


8000 


2.2 


10.5 


4.8 


NGC 4388 


0.0084 


170 


2001-06-08 


A-S 


20 


6332 


W2/FR 


280 


3.4 


13.4 


3.9 


NGC 4507 


0.0118 


240 


2001-03-15 


A-S/H 


140 


8259 


W2/FR 


520 


3.3 


11.0 


3.3 


NGC 5347 


0.0078 


160 


2004-06-05 


A-S 


37 


8598 


W2/FR 


800 


0.27 


0.8 


3.0 


NGC 5643 


0.0040 


80 


2004-12-26 


A-S 


8 


5411 


W2/F5 


700 


1.4 


5.1 


3.6 


NGC 7212 


0.0266 


530 


2003-07-22 


A-S 


20 


6332 


W2/FR 


600 


0.8 


8.8 


11.0 



" Radial extent of 1 arcsec at the distance of the source, assuming the cosmology reported in the text 
* Chandra instrument (A-S: ACIS-S; A-S/H: ACIS-S/HETG) 

Chandra exposure time in ks 
'I HST instrument/filter configuration (W: WFPC; W2: WFPC2; F: FOC; F5: F502N; FR: FR533N) 

HST exposure time in s 
^ All values taken from this paper, except for NGC 4388 ilwasawa et all2003l) 
^ All values taken from lSchmitt et alji20033) . except for NGC 5643 (this paper) 



servations of Mrk 3, using the standard SAS ( iGabriel et alJ 
[2004) meta-task rgsproc, and the most updated calibration 
files available at the moment the reduction was performed (July 
2005). Each spectrum was produced assuming that the X-ray 
emission comes from a point source located at the optical co- 
ordinated of the Mrk 3 AGN, thus ensuring an homogeneous 
wavelength reconstruction. The wavelength systematic uncer- 
tainty is 8 mA across the whole RGS sensitive bandpass. 
Spectral fits have been performed on two 1 st order RGS 1 and 
RGS2 spectra, resulting from the sum of all the available spec- 
tra, generated by rgscombine. Fig.0]- which we present for 
illustrative purpose only in Sect. 13.2.21 - was generated by com- 
bining all the 1st and 2nd order spectra of the two RGS cam- 
eras, and smoothing the combined spectra through a convolu- 
tion with a 5-spectral channel wide triangular function. 

A ll spectra were analyzed with XSPEC 12.2.0 jArnaudl 
Il996l) . In the following, errors correspond to the 90% confi- 
dence level for one interesting parameter (A^^ = 2.71), where 
not otherwise stated. 

3. Analysis 

3.1. Imaging 

In none of the sources there is evidence of hard X-ray emission 
extended beyond the instrument Point Spread Function (PSF), 
in agreement with an origin of this component within a few 
pc from the nucleus. On the other hand, most of the objects 
present a clear extension in the soft X-ray emission, with 50 to 
80% of soft X-ray counts detected outside the inner 1 arcsec 
(which corresponds to 60-530 pc, depending on the source: see 
third column of Table 0. These percentages are much higher 
than those expected for an on-axis point-like source observed 
by ACIS (see Fig. 4.7 and Table 4.2 of The Chandra Proposers' 
Observatory Guide'). The soft X-ray emission is much more 

' http://cxc.harvard.edu/proposer/POG/html/MPOG.html 



concentrated in NGC 7212, where only 25% of the counts are 
detected outside the inner 1 arcsec. This is not surprising, as 
this source is by far the farthest of our sample (see Table 
However, the emission is strongly elongated and so clearly ex- 
tended. Finally, Matt et al. (2004) report no extension in the 
Chandra image of NGC 4507. We also do not find any con- 
clusive evidence in favour of an extended soft X-ray emission 
in our re-analysis of the data. However, it should be noted that 
the observation is severely affected by pile-up (see Sect. 12. 2> . 
which prevents us from performing a quantitative comparison 
with the PSF of the instrument. 

Figure ^ shows the contours of the Chandra soft X-ray 
emission superimposed on the HST [O III] images, for all the 
sources in our sample. Since we are interested in the relation- 
ship between the X-ray and the [O III] morphologies, we tried 
to adopt a similar dynamical range to display the two compo- 
nents, based on the peak flux of the image. The choice of the 
dynamical range to use for each source was basically dictated 
by the brightness of the nucleus in comparison to the weaker 
extended emission. This criterion, though being just one of the 
possible choices, allows us to have fairly homogeneous images 
for all the sample. 

The coincidence between the soft X-ray and [O III] emis- 
sion is striking, both in the extension and in the overall 
morphology. Unfortunately, the lower spatial resolution of 
Chandra with respect to HST does not allow us to perform a 
detailed comparison of the substructures apparent in the latter. 
Moreover, the uncertainty in the relative alignment of the opti- 
cal and X-ray images (see Sect. 12.2.21 contribute to hinder this 
comparison, which would indeed be crucial to understand the 
real nature of the two components (see Sect. |5}. Two sources 
deserve some further comments. In the case of NGC 4507, as 
already mentioned, the soft X-ray extension is not clear: how- 
ever, the [O III] emission is also extended on a very small scale 
and rather symmetrically. As for NGC 4388, it should be noted 
that the X-ray emission (which extends at least up to 16 kpc 
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Fig. 1. Chandra soft X-ray contours superimposed on HST [O III] images. The contours correspond to five logarithmic intervals 
in the range of 1.5-50% (NGC 1386), 5-80% (Mrk 3), 5-90% (NGC 3393), 4-50% (NGC 4388), 0.5-50% (NGC 4507) and 
2-50% (NGC 5347) of the peak flux. The HST images are scaled with the same criterion for each source, with the exception of 
NGC 4388 and NGC 4507, whose [Olll] emission goes down to the 2% and the 0.1% of the peak, respectively. The images are 
aligned as explained in Sect. I2.2.:?l and centered on the Chandra count peak. 
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Fig. 2. Same of Fig.[T] but for NGC 5643 (8-80%) and NGC 7212 (1-50%). In the case of NGC 5643, the [O III] emission goes 
down to the 0.5% of the peak. 



at the lowest energies) may well be influenced by the Virgo 
cluster where the AGN lies. Moreover, this object may also 
be peculiar with respect to the others in our sample, because it 
presents also a Fe Ka extended emission ( I wasawa et al..2003) . 

Generally, the observed structures extend up to scales of 
kpc. While the morphologies for most sources could be ap- 
proximately described as a bi-conical (or at least mono-conical) 
structure, for some objects a more complex representation is 
needed. In particular, Mrk 3 and NGC 3393 present S-shaped 
arms, which cannot be understood simply as collimation of the 
nuclear radiation by a torus. We will discuss the role of radio 
jets in the formation of these substructures in Sect.|5l 

3.2. Spectral analysis 

As already mentioned in Sect. 12.2.11 we analyze in this pa- 
per for the first time the ACIS-S spectra of 5 out of the 8 ob- 
jects in our sample. For the other three sources (namely Mrk 3, 
NGC 4388 and NGC 4507) we will refer to the results on their 
Chandra spectra which have been already published in litera- 
ture. However, in the case of Mrk 3, we will also present new 
results from a combined 190 ks XMM-NewtonfRGS spectrum. 

In all the sources where the statistics were good enough to 
perform a comparison between the soft X-ray spectra from the 
nuclear and the extended region (as defined in Sect. 12.2. 1> . no 
significant differences were found between the two. In partic- 
ular, the observed emission Unes are required in both spectra. 
Therefore, all the following fits were performed on the total 
spectrum, which includes counts from both regions. 

Two different models were tried to fit the soft spectra an- 
alyzed in this paper: a 'thermal' model, reproducing emis- 
sion from one or more optically thin, collisionally ionized 
plasma with free elemental abundances (Mekal in XSPEC: 

I'- — -» F 
Mewe et al. 1985); a 'scattering' model, i.e. a power-law and as 

many unresolved emission lines as required by the data, repro- 
ducing reflection of the nuclear continuum from a Compton- 



Table 2. Fluxes and iron line properties for the Chandra ob- 
servations analysed in this paper. The iron line fluxes and EWs 
refer to the 'local' fits (see text for details). 



Source 


Fo.5- 


a 

2keV 


F2-10kcV " 


EWpe* 


Ipe 


C 


NGC 1386 


l.f 


'-0.5 


2.4 ± 1.0 


1.5 


+0.8 
-0.5 


6+"^ 

"•"-0.2 


NGC 3393 


2.: 


)+0.5 
--0.2 


2.6 ± 0.8 


0.6 H 


= 0.3 


0.4 ± 


0.2 


NGC 5347 


0.27 


±0.06 


2.4 ± 0.7 


1.7 H 


= 0.7 


0.5 ± 


0.2 


NGC 5643 


1.4 


±0.2 


6 3+2' 

U.J 2.0 


1.5 


+0.8 
-0.6 




NGC 7212 


O.f 


;+0.1 
-0.2 


5 2+'-^^ 

-1.2 


0.7 H 


= 0.3 


0.7 ± 


0.3 



"In units of 10"'^^ erg cm"^ s"' 
*Iron line EW in units of keV 
''Iron line flux in units of 10"^ ph cm"- s"' 



thin material photoionized by the AGN itself. In the latter 
model, the photon index of the powerlaw is left frozen to the 
value of r = 1 .7, which is a typical value for the primary con- 
tinuum of Seyfert galaxies. 

In the hard band, all the sources analysed in this paper are 
dominated by the iron emission line, while statistics are in most 
cases too poor to give significant constraints on the underly- 
ing continuum. However, when an additional component is re- 
quired to model the hard X-ray spectrum, it will be explicitly 
reported in the text. As a consequence of the quality of the hard 
X-ray spectra of the sources, the iron line properties have been 
measured through the 'local' fits (described in Sect. |^}: the 
values reported in Table|2]refer to these fits. 

Detailed results on each source are described in the follow- 
ing subsections. 
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Table 3. NGC 1386: emission lines included in the fit of the 
soft X-ray spectrum, along with their most likely identifications 
(see text for details). 



Energy 

^Kev J 


Flux 

\ L\J pil Clll S ) 


Id. 


53+° '" 


^■"-1.5 


N VII Kq' - O VII Kq" 


U.OD ± U.UZ 


1 Q 1 n A 
1 ± u.^ 


KJ Vlll JS.U^ - Vii J\p 


0.76 + 0.01 


1.5 ±0.4 


VII RRC - VIII K/3 




1.5 ±0.3 


OVIIIRRC 


94+0 '" 


1.1 ±0.3 


Ne IX Ka 


1.05+0.02 


0.6 ±0.2 


Nex Ka 


1 16+°''^ 


0.5 ±0.2 


Fe XXIV L 


1.38 + 0.04 


0.2 ±0.1 


Mg XI Ka 


1.78 ±0.05 


0.2 ±0.1 


Si XIII Ka 



3.2.1. NGC 1386 

The source was found to be a reflection-dominated object in 
the XMM-Newton observation ( Guainazzi et al. 2005a.) . This 
is confirmed by the strong iron line in the Chandra spectrum 
(see Fig. |3] and Table |3, which by far dominates the hard 
spectrum and whose properties are fully consistent with those 
found with XMM-Newton. Moreover, both the soft and the 
hard X-ray flu xes are in agreement with the ones reported by 
lOuainazzi et al.. (.2005a) . 

The EPIC soft X -ray spectra were best fitted with two 
thermal components fcuainazzi et alJ Eo05al) . However, this 
model, though leading to values similar to those found with 
XMM-Newton (temperatures of 0.14+" "5 ^nd 0.72+°{;^ keV, 
with a metal underabundance of O.O?;^!] !]^), fails to produce a 
good fit on the Chandra data (x^ = 62/38 d.o.f.). A much bet- 
ter fit is achieved with the 'scattering' model (x^ = 18/20). The 
lines included in the model are from H- and H-like N, O, Ne, 
Mg, possibly blended with Radiative Recombination Continua 
(RRC) from C and O (see Table|3}- The latter features, if con- 
firmed by future, high resolution spectra of NGC 1386, would 
strongly favour an origin of this emission from photoionized 
plasma, like in the cases of NGC 1068 ( K inkha^ walag^D 
2002; Brinkman et al.' l2002l) and Mrk 3 llSako et alJ 1200(1 
Bianchi et al. 2005b) . 

A final comment must be spent on a feature found at 0.29 
keV. While it is possibly associated with the C V Ka triplet, its 
large flux (of the order of 10"^ ph cm"^ s"') makes its identi- 
fication as an emission line troublesome. However, the feature 
lies at the end of the band considered 'good' for ACIS spec- 
tra: below 0.3 keV calibration and background issues prevent 
us from performing a reliable analysis. 



3.2.2. Mrk 3 

The Chandra ACIS-HETG and the XMM -Newton EP I C and 
RG S spectra wer e analyz ed in detail by ISako et all tOOd) 
and ' Bianchi et all (|2005b'), respectively. It is one of the few 
Compton-thick Seyfert galaxies that have high resolution spec- 
tra, showing that the soft X-ray spectrum is indeed dominated 
by emission lines. The ratio between the forbidden to the reso- 
nant line in the He-like triplets of light metals and the presence 
of RRC from the same elements give strong evidence in favour 
of an emitting gas photoionized by the nuclear continuum. 

We show in Fig. I^the high resolution spectrum with the 
largest S/N available for Mrk 3 and one of the best high res- 
olution spectra for an AGN ever measured. It was obtained 
summing all the RGS observations of the source (for a total 
of ^ 190 ks) and combining RGS 1 and RGS 2 first and second 
orders. 

The spectrum is dominated by K-shell emission lines from 
O and Ne, with weaker Fe L lines mainly from Fe XVII and 
Fe XVIII. Spectra from plasma in collisional ionization equi- 
librium (CIE) are instead dominated by Fe L emission line s 
(see e.g. the spectrum of solare flares: fPhilUps et alJll982i) . 
As already observed bv lSako et al.l ll2000l) . this could in prin- 
ciple be accounted for by an at least one order-of-magnitude 
iron underabundance, but this is at odds with the EW of the 
iron Ka line and the depth of the neutral iron edge, which are 
consistent with an iron a bundance only slightly lower than the 
solar one (Bia nchi et al ]|2005b). Moreover, the clear detection 
of narrow O VII and O VIII RRCs implies an electron temper- 
ature much cooler than the o ne expected in a CIE plasma with 
the sa me ionic species (e.g. iLiedahl & Paerelslll996t iLiedahll 
Il999l) . Therefore, the gas is likely to be in photoionization equi- 
librium (PIE). 

Indeed, the spectrum is telling us something more. The 
resonance lines in the O VII and Ne IX triplets are stronger 
than those expected fo r pure recombination in PIE (e.g. 
IPorauet & Dubaul l2000l) . The resonance 3d-2p Fe L lines 
should also be suppressed when recombination is the domi- 
nant process (Liedahlet al. 1990), while these transitions are 
clearly detected in the observed spectrum together with 3i- 
2p lines. Both problems can be solved taking into account 
the role of reso nant scattering (.Band et al. 1990: Matt 1994; 
iKrolik & Krisslfl995i) . Resonant lines in He-hke triplets are 
greatly enhanced by this mechanism with respect to recombina- 
tion, es pecially for low column densities (see e.g. Bianchi et al. 
l2005ah . Similarly, photoexc itation significantly f avors the pro- 
duction of Fe L 3d-2p lines ('Kall man et al.ll996l) . 

While a complete and detailed analysis of the combined 
RGS spectrum is beyond the scope of this paper, we com- 
pared the O VII-O VIII series with the photoionization code 
Photoion ( Kinkhabwala et al. 2003) , in order to get some 
more quantitative information from this high S/N part of the 
spectrum^. We first only dealt with the spectrum between 0.55- 

^ As already reported by iKinkhabwala et alj 12002*), PHOTOION, 
though being a local model in XSPEC, does not allow a real fit proce- 
dure, because its complex nature makes it very difficult to deal prop- 
erly with error calculations. Therefore, our 'best fit values' have been 
found by trial and error, and we cannot estimate an error 
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0.66 keV (observer's frame), in order to include only the O VII 
triplet and the O VIII Lya. The strong resonant line with re- 
spect to the forbidden and the intercombination lines is best 
reproduced with a O VII column density of 7 x 10'^ cm"^ and a 
transverse velocity distribution o""'' = 250 km s"' . Indeed, both 
parameters contribute to make resonant scattering an impor- 
taiit process in the gas under analysis (e.g. Kinkhabwala et al. 
l2002t feianchi et aimOOSal) . The observed forbidden line ap- 
pears stronger than expected with respect to the intercombina- 
tion line: this was also observed in the spectrum of NGC 1068 
iKinkhabwala et al . 2002; B rinkman et alJ2002l) and explained 
with an additional contribution due to inner-shell photoioniza- 
tion in O VI (see also Kinkhabwala et al. 2003, for further de- 
tails). The strong O VIII Lya line requires a much higher col- 
umn density for this ion (1.6 x 10"^ cm"^), thus suggesting that 
the average ionization parameter could be high. 

The inclusion of the spectrum up to 0.79 keV (observer's 
frame) shows that the flux of lower order lines (O VII and O VIII 
yS) and the O VII RRC is underpredicted by the parameters used 
to account for the a lines. Again, this is another piece of evi- 
dence in favour of a strong contribution from photoexcitation 
and, in general, an important in dicator of photoionization as 
the dominant process at work ("see lKinkhabwala et alJ200^ for 
a detailed discussion on this issue). In order to reproduce cor- 
rectly also these lines, we need lower column densities (3x10'^ 
and 8 x 10'^ cm"^ for Ovil and Ovill, respectively) and a 
lower o-\^'^ = 200 km s"' . Note that it is not easy to disentangle 
the role of each of the two parameters separately, since both 
play on the relative role between recombination and photoexci- 
tation. With these parameters, the need for an additional contri- 
bution for the O VII forbidden line flux is higher Moreover, the 
observed O VII RRC keeps on being stronger than the expected 
one. We do not have an easy explanation for this behaviour, but 
it is likely that the close Fe XVII 3s-2p line (which is not mod- 
elled by Photoion) may contaminate the flux of the RRC. 

In summ ary, we c onfirm previous results by ISako et all 
(I2OOO) and Bianchi et al. ( ,2005b,) . but on stronger evidence: 
the soft X-ray spectrum of Mrk 3 is produced in a gas in pho- 
toionization equilibrium, with a significant contribution from 
resonant scattering as well as recombination. 

Moreover, visual inspection of Fig.|4]also suggests that the 
gas is not strongly in/outflowing, since large shifts of the emis- 
sion lines are not observed with respect to the theoretical val- 
ues. This is quantitatively supported by the results on the spec- 
tral fits, which lead to tight constraints on the shifts of the O VII 
forbidden line, the O VIII Ko' and the N VII Ka, being respec- 
tively + 60, 30 + 1 10 and 10+}^" km s"' (note that the system- 
atic uncertainty of the RGS at these energies is around 100-130 
km s"'). However, this may be simply due to the fact that we 
are observing a spectrum which includes both sides of the bi- 
cone, thus naturally producing a zero net shift. In any case, the 
line widths are also all upper limits, being < 300, < 400 and 
< 1 000 km s ' for the above-mentioned lines, respectively. 



Table 4. NGC 3393: emission fines included in the fit of the 
soft X-ray spectrum, along with their most likely identifications 
(see text for details). 



Energy 

(keY) 


Flux 

(10^^ ph cm"- s"') 


Id. 


041+0.04 

"■^^-0.06 


2 8+2-^ 

Z.O_23 


r V RRr - N VT Kry 


"■^°-o.oi 


4.6 ± 1.2 


N VTT Kry 

ly V 1 ± ±^t_i: 


0.573 ± 0.001 




VII Ka 


"■"'-0.01 


9 9+0.6 
^•^-0.5 


VIII Ko- - VII igS 


77+0 0' 


2 5+°-^ 
^•-'-0.4 


VII RRC - VIII K/3 


89+""" 

"■°'^-0.02 


3.0 ±0.4 


VIII RRC - Ne IX Ka 


1 03+"°' 
^■"-^-0.02 


9+" ' 


Nex Kff 


1 16+" °2 


3+° ' 

"•-'-0.2 


Fe XXIV L 


1 30+°°2 


0.3 ±0.1 


Mg XI Ka 



3.2.3. NGC 3393 

On the basis of its BeppoSAX and XMM-Newton observa- 
tions, G uainaz zi et al. ( 2005a ) concluded that NGC 3393 is 
a Compton-thick source, even if the iron line was barely de- 
tectable because of low statistics. Our local fit on the Chandra 
spectrum confirms the presence of a strong iron line as the dom- 
inant component of the hard spectrum (with a flux consistent 
with the one observed in the EPIC spectra), in agreement with 
the reprocessing-dominated nature of the object. 

We first tried to fit the soft X-ray spectrum with a 
thermal model, as suggested by the XMM-Newton analysis 
JCiuainazzi et al. 2005a). We get parameters very similar to 
those found in the EPIC spectra (two Mekal components, 
with temperatures of 0.17;^qq2 and 0.70;^qq3 keV and an un- 
derabundance of 0.06 ^^JJ^j with respect to the solar one), but 
the fit is unacceptable (x^ - 103/58 d.o.f.), mainly because 
of large residuals at the O VII and Ne IX Ka energies. We get 
a much better fit with a photoionized model (x^ - 33/44). 
The required lines are similar to the ones found for NGC 1386 
(see Table 0. The presence of a line at ^ 1.8 keV, detected 
by Guainazzi et al. (2005a) in the XMM-Newton spectrum and 
expected in our model as emission from Si XIII Ka, is not re- 
quired by the Chandra data, being a detection lower than 2 cr. 

3.2.4. NGC 4388 

llwasawa et alJ (l2003l) analyzed in detail both the image and 
the spectrum of the Chandra observation of this source. They 
found a strongly absorbed (Nh - 4 x lO^-' cm"^) spectrum, 
together with an extended soft X-ray emission, whose most 
likely origin is in a gas photoionized by the central AGN. 
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Interestinglv. llwasawa et al suggested a common origin 

for the [O III] and the soft X-ray emission. Moreover, they also 
reported that no significant variations are found in the spectral 
shape along the radius of the extended emission, so that the 
ionization parameter remains fairly constant, thus implying a 
density which decreases with the distance roughly like r"^. 

3.2.5. NGC4507 

The Chandra spectrum of NGC 4507 was analyzed in detail by 
iMatt et al . (2004), together with the XMM-Newton data. The 
source is seen through a large column density {Nh - 4 x 10^^ 
cm"^), while the presence of a Compton-thick reflector is in- 
ferred by the presence of a Compton reflection component, a 
large iron Ka line and its Compton shoulder. The soft X-ray 
emission is dominated by emission lines, likely produced in 
a photoionized material. Interestingly, Matt et al. (2004) found 
also an absorption line from Fe XXV, thus suggesting also the 
presence of an highly ionized gas, as observed in a large num- 
ber of Seyfert galaxies (see e.g. .Bianchi et aL.2005at and ref- 
erences therein). 

3.2.6. NGC 5347 

The Chandra spectrum of this object is clearly reflection- 
dominated: a good fit (y^=9/12 d.o.f.) is achieved with a 
Compton reflection component, a steep (F ^ 2.4) powerlaw 
and a neutral iron line with a large EW (see Table |2}. The low 
statistics do not allow us to analyze in detail the soft X-ray 
spectrum, which can be also modelled by two thermal compo- 
nents, with loosely constrained temperatures around 0.16 and 2 
keV and quasi-zero abundances {x^=%l9 d.o.f.). 

3.2.7. NGC 5643 

The XMM-Newton observation suggested a complicated pic- 
ture for NGC 5643, with a heavily absorbing column den- 
sity which likely obscures part of its own Compton reflec- 
tion at high energies and a soft X-ray spectrum dominated 
by mainly unresolved emission lines produced in a photoion- 
ized plasma. Other pos sibilities cannot be, however, ruled out 
(^Guaina zzi et al J 12004). The hard X-ray flux measured with 
Chandra is consistent with the one found with XMM-Newton, 
while the soft fluxes, though marginally, are not in agreement. 

The quality of the Chandra spectrum of NGC 5643 is 
not good enough to perform a detailed analysis in the soft 
band. However, a reasonable broad band fit (x^=5l\0 d.o.f.) is 
achieved with a pure Compton reflection component, an iron 
line and a steep (F > 3) powerlaw to model the softer part of 
the spectrum. The strong iron line is confirmed by the local fit 
(see Table |2j- In this source, the statistics do not allow us to 
check if this component may be reproduced by the blending of 
the emission lines expected from a photoionized plasma, even 
if the 'KMbA-Newton an alysis revealed that this is indeed pos- 
sible JGuainazzi et alJ2004.) . However, the thermal model with 
two Mekal components gives a fit of similar quahty {x^='in 



d.o.f.), with temperatures of 0.15 and 0.62 keV and a metal 
abundance lower than 0.08 times the solar. 

3.2.8. NGC 721 2 

The soft spectrum is well fitted by the scattering model with 
two emission fines (at 0.86+|5 |55 and 1.84+° °^ keV) required by 
the data at the 93% and 95% confidence level according to F- 
test, likely identified as Ka emission from Ne IX and Si XIII. 
Indeed, the bad modelling of these residuals is the main reason 
why a thermal model fails to provide a good fit (x^ - 18/13 
d.o.f.). The broadband fit includes also a pure Compton reflec- 
tion component and a strong iron line, confirmed by the local 
fit (see Table |2ll, thus suggesting that the source is Compton- 
thick. This is in agreement with the XMM-Newton observation 
(Guainazzi et al. 2005b), whose soft and hard X-ray fluxes are 
also consistent with those measured by Chandra. 

4. Photoionization models 

The spectral analysis of the sources in our sample suggests that 
their soft X-ray spectra are likely dominated by emission lines 
produced in a material photoionized by the central AGN. On 
the other hand, the striking resemblance of [Olll] structures 
with the soft X-ray emission favors a common origin for both 
components. Therefore, since the NLR is generally believed to 
be also produced mainly by photoionization, we generated a 
number of models in order to investigate if a solution in terms 
of a single photoionized material to produce the optical NLR 
and the soft X-ray emission is tenable. 

Calculations were pe rformed with version 96.01 of 
Cloudy, last described bv lFerland et alJ (Il998h . The adopted 
model is represented by a conical geometry, which extends 
from 1 to 350 pc from the nucleus, with temperature set by 
photoionization equ ilibrium under a typical AGN continuum 
jKorista et alJll997i) . We produced a detailed grid of models, 
as a function of the following parameters: ionization parameter 
U of the illuminated face of the gas, defined as: 



U 



Jvb hv 



4;rr2, 



(1) 



where c is the speed of light, r the distance of the gas from 
the illuminating source, n^, the electron density and v« the fre- 
quency corresponding to 1 Rydberg; filling factor of the gas 
/; density of the gas «f . We assume a power-law radial depen- 
dence of the density, njro) ^, with p - Q (constant den- 
sity) and varying between 1 and 2.4. We limit our models to a 
minimum total column density of 10^" cm^^. Note that the ra- 
dial column densities under analysi s are always consistent \y ith 
a negligible edge opacity (see e.g. Kinkhabwala et al.l20o3) . 

The ratio between the [O III] /15007 line and the soft X-ray 
emission (defined as the total flux of the Ka and KjS emission 
lines from N, O, Mg, S, Si in the range 0.5-2.0 keV) was cal- 
culated for each set of parameters. In Fig.|5] each symbol rep- 
resents a solution in the U versus plane, where this ratio has 
a value within 2.8 and 11, as observed in our sample (see Table 
[0. Different symbols correspond to different values of yS. The 
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Fig. 5. Each symbol represents one solution in our grid of 
Cloudy models that satisfies the condition of total [O iii] to 
soft X-ray flux ratio in the range observed in our sample, plot- 
ted in a three-parameter space U, rig and /3, i.e. the ionization 
parameter and the density at the beginning of the cone of gas 
(1 pc), and the index of the density powerlaw, represented by 
different symbols {triangles: fi — 1.6, diamonds: yS = 1.8, stars: 
j3 = 2.0, squares: /3 = 2.4). The horizontal lines determine the 
limit corresponding to a total column density of 10^** cm"^ for 
each index: solutions below this limit are not plotted. Solutions 
with y6 - 2.2 are not plotted for clarity reasons. The net effect 
of the filling factor / is to shift the density values: the two y 
axes refer to / = 0.1 and / = 0.01 (see text for details). 

net effect of changing the filling factor is simply a shift of the 
solutions along the density axis, by a factor equivalent to the 
variation in /, thus reproducing the same total column density 
for each set of three parameters consituting a 'good' solution. 
The solutions occupy well-defined regions in the ng - U dia- 
gram, with those with lower /3 being at larger values of ioniza- 
tion parameters. 

The reason for this behaviour becomes clear inspecting Fig. 
|6l where the [O III] to soft X-ray flux ratio is plotted as a func- 
tion of the radius of the gas. Since U oc n^^r^^, all density 
laws with /? < 2 produce a gas with a ionization parameter de- 
creasing along with the distance. In these cases (/? = 1.6 and 
1.8 in Fig.|6}, most of the soft X-ray emission is produced in 
the inner radii of the cone, while the bulk of the [O III] emis- 
sion is produced farther away, where the gas is less ionized. If 
13-2, the ionization parameter remains fairly uniform up to 
large radii, so that the total observed ratio between [O III] and 
soft X-ray is constant with radius. Finally, if j6 > 2, the ioniza- 
tion parameter increases with radius, so most of the soft X-rays 
are actually produced at larger radii, while the [O III] emission 
line is mostly concentrated around the nucleus (cases [3 - 2.2 
and 2.4 in Fig.|6j. 

This is the reason why solutions with < 1.6 and (3 > 
2.4 are not plotted in Fig.|5] even if these exist, satisfying the 
overall [O III] to soft X-ray flux ratio. Their radial behaviour 



Fig. 6. The [O III] to soft X-ray ratio plotted as a function of the 
radius of the gas, for different values of /3, when log ne - 3 and 
/ = 0.1. The corresponding values of log U for these solutions 
are (from top to bottom): -0.85, -1.15, -1.4, -1.5, -1.6. 



is radically different from what seen in Fig. QandlJ] where 
both the [O III] and the soft X-ray emission are produced up to 
large radii. In particular, it is worth noting that constant-density 
solutions are totally unacceptable. 

Before going through the discussion on the results of these 
calculations, we would like to point out some caveats. The flux 
included in the quantity we call 'soft X-ray' does not include 
Fe L emission, which is known to be quite important, even if 
not dominant, in photoionized spectra (see, for example. Fig. 
13. Moreover, no contribution from the scattered continuum has 
been taken into account in the total soft X -ray flux. However, 
this should not be very re levant (see e.g. iKinkhabwala et all 
'2002*: 'Bia nchi et alJ2005 bV In any case, the inclusion of these 
components is not expected to change dramatically the value of 
the [O III] to soft X-ray ratio. 

We have not investigated the effects of varying the incident 
continuum. In particular, some authors include X-ray and UV 
absorbers between the nucleus and the NLR clouds, finding ev- 
idence that an absorbe d incident contin uum is actuall y needed 
by th e data (see e.g. [Alexander et aL 1999; Kraemer et alJ 
hood) . The presence of this absorber would change the number 
of ionizing photons in the incident radiation on the NLR clouds, 
thus effectively changing the ionization equilibria calculated in 
our model. We will discuss the role and the properties of these 
absorbers in the next section, particularly regarding their pos- 
sible link with the gas analysed in this work. 

Finall y, the role of dust wa s not included in our models. 
However, Galliano et alJ ( l2003l) found that the mid-infrared 
(MIR) emission in NGC 1068 is very well correlated with the 
[O III] NLR, thus suggesting the presence of dust grains in this 
material. The dust could have relevant effects on the absorp- 
tion of optical emission, thus changing the value of the ratio 
between [O III] and the soft X-ray flux. 
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As a consequence of these caveats, these models should be 
treated as simple toy models. In order for more complex models 
to be tested, one would need X-ray imaging detectors match- 
ing the HST sub-arcsec resolution and possibly space resolved 
high-resolution X-ray spectroscopy at the same scale, but both 
requirements will not be available for many decades to come. 
Our results are a proof of existence of simple scenarios consis- 
tent with the currently available experimental data. However, 
they do not imply that more complex or different scenarios can 
be ruled out, as we will discuss in detail in the next Section. 

5. Discussion 

The analysis of the Chandra images presented in this paper 
has unambiguously shown that there is a spatial correlation 
between the soft X-ray emission in Seyfert galaxies and their 
[Olll] emission. Moreover, the spectral analysis of the same 
sources suggests that the origin of the soft X-ray component 
is likely in a photoionized gas. In one case, Mrk 3 (but see 
also the cases of NGC 1068 and Circinus, not included in this 
sample), this component is actually resolved in plenty of emis- 
sion lines in high resolution spectra and a number of diagnostic 
tools indicate that they are produced by gas in photoionization 
equilibrium. In all the other sources, where only CCD spectra 
are available, a thermal model is always disfavored, in most 
cases, on statistical grounds, but also because it leads to un- 
physically low metal abundances. In the latter case, however, 
a thermal component cannot be completely ruled out, since 
low abundances can be, in principle, mimicked by iron deple- 
tion from dust in interstellar clouds ( Ptak et al..i.l99 7) or dilu- 
tion by featureless con tinua, produced, for example, in binaries 
jGuainazzi et alJl999l) . 

The photoionization models presented in the previous sec- 
tion suggest that it is in principle possible to have a gas, pho- 
toionized by an AGN and extending for hundreds of pc, which 
emits the right ratio of [Olll] to soft X-ray with a reason- 
able radial profile. These results imply peak densities around 
10^ cm"^ (for a filling factor of 0.1), which are reasonably 
in agree ment with wh at found for the NLR of many sources 
(see e g.lCooke et al.lEo OO: Nelson et al. 2000; FraauelH et al] 
l2003t iBradlev et alJl20041) . Moreover, solutions with constant 
density must be rejected, and density laws that fall similarly to 
r"^ have to be preferred. 

We plot in Fig. 0the ratio between the HST and the soft 
X-ray Chandra counts as a function of the distance from the 
nucleus for the source with the best Chandra data in our sam- 
ple, NGC 3393. The counts at each radius are taken from rect- 
angular regions 0.5 arcsec wide and 5 arcsec long, perpendic- 
ular to the axis of the extended region. The HST image was 
convolved with a Gaussian with FWHM=0.8 arcsec to mimic 
Chandra resolution. The plot shows that this ratio is fairly 
constant within a factor of a few, thus requiring an ioniza- 
tion parameter which changes slowly with radius. In partic- 
ular, the data points in Fig. suggest a decrease somewhat 
flatter than 2, but still steeper than r ' Interestingly, sim- 
ilar results are found for the NLR on the basis of the ratio 
of the [Sill] ^6717 to ^6731, [Olll] ^5007 to H£and_[Olll 
/13727 to HyS for a number of sources, (see e.g. lKraemer et alJ 
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Fig. 7. NGC 3393: HST to Chandra soft X-ray counts ratio as a 
function of distance from the nucleus. The two curves refers to 
the photoionized models plotted in Fig.|6l(see text for details). 

bOOOHKraemer & Cre nshaw'200d'Bra dlev et alJ2004l) . includ- 
ing, notably, Mrk 3 (Collins et al. 2005). 

A simple explanation for this density decrease with radius 
may be that the NLR clouds are not pressure confined an d so 
they tend to expand outwards, as suggested by Kraemer et alJ 
(2000). A density fall like r^^ would arise naturally if the 
outflow had a constant velocity and mass rate. However, 
this picture is clearly oversimplified and other mechanisms 
can contribute to the observed radial behaviour of the ion- 
ization parameter. In particular, the [Olll] structures are 
generally found very well correlated with the radio images, 
thus suggesting an important role of jets in the overall NLR 
shape. We refer the reader to the papers where the radio 
images of these sources are presented to see the correlation 
with the [O lllVsoft X-rays component s shown in Fig, [p 
NGC 1386 llNagaretalJfl99"9i) Mrk 3 llCapetti et alJll99i)! 
NGC 3393 llCooke et alJ l2000l). NG C 4388 ( iKukula et alJ 
1995* ), NGC 4507 (' MorgantietalJ Il999l). NGC 5643 
( Morris etal J ^985) . NGC 7212 JPalcke et alJ Il998l) . On 
the other hand, it is worth noting that NGC 5347 represents an 
exception, since it does not pr esent extended radio emission 
jThean et alJ200(ilSchmitt et al.i2003b.) . 

In particular, it has been proposed that the expansion of the 
radio ejecta compresses the gas, enhancing its radiative emis- 
sion, thus effectively dominating the observed morphology of 
the NLR. Moreover, Sutherland et al. (1993) have also sug- 
gested that the interaction between the radio plasma and the 
surrounding medium produce an ionizing continuum capable 
of photoionizing the same gas. In other words, the local ion- 
ization parameter in a particular region of the NLR may be 
affected by this interaction as well as from the AGN contin- 
uum. It is clear that this phenomenon would reduce the need 
for a steep decrease of the density with the radius, since local 
sources of photoionization would prevent the ionization param- 
eter from falling steeply with radius. 
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We have already mentioned the presence of the X-ray ab- 
sorber in the previous section. The so-called 'warm absorbers' 
are often observed in the spectra of Seyfert Is and are thought 
to be located between the Broad L ine Region (BLR) and the 
torus (see e.g. iBlustin et al.ll2005l and references therein). If 
this is the case, this material is not supposed to be an important 
emitter in Seyfert 2s or not even be observable because ob- 
scured by the torus. On the other hand, the gas under analysis 
in this paper may well be along the line of sight when observing 
a Seyfert 1 and so effectively add its column density to the in- 
ner warm absorber. The column densities of warm absorber are 
measured approximately in the range 10^' - lO^"* cm"^, but are 
believed to be composed of various phases, possibly with dif- 
ferent geometrical and kinematic properties (e.g. lBlustin et alJ 
l2005h . In this sense, it is still possible that the denser and inner 
part of the soft X-ray emitting region resulting in the scenario 
discussed in this paper may represent the outer part of the warm 
absorber. 

In this regard, it is interesting to compare the gas velocity 
as inferred from the RGS spectra of Mrk 3 with those typically 
found in the warm absorbers and in the NLR. As we have seen 
in Sect. 13.2.21 the energy shifts measured in the X-ray emission 
lines of Mrk 3 imply that the gas has a very low bulk veloc- 
ity, no larger than around 100 km s On the other hand, the 
warm absorbers are often found to be ou tflowing at velocitie s 
larger than a few hundreds km s ' (e.g. iBlustin et alJ l2005). 
This is, at least qualitatively, in agreement with the velocity 
maps measured for the NLR, which resemble those of the warm 
absorber in the region closer to the nu cleus, then decreasing 
toward zero velocity outwards (see e.g. iRuiz et al ]|2005 t and 
references therein). 

Indeed, it has been suggested that even the NLR could be 
made up of different components, with different density and 
ionization (e.g. Korista & Ferland 1989; Kraemer et al. 2000). 
We have clear evidence of the need of materials with differ- 
ent ionization stages also in X-rays, with the observation, in a 
large number of objects, of lines from highly ionized iron, in- 
consistent with the other line s from lighter metals observed in 
the soft X-ray spectra (see e.g. lBianchi et alJ2005 al). Moreover, 
all the solutions presented in this paper have negligible emis- 
sion from Si, while H and He-like emission from this element 
are detected even in some of the sources included in our sam- 
ple. In this sense, it is likely not surprising that the values of 
U of our solutions are larger than the or ies tvpically required 
by t he modelling of the NLR (see e.g. iKraemer et al.lll999[ 
iBas kin & Laor 2005). The latter are probing the gas dominated 
by the optical emission, while ours are basically driven by the 
need of a strong X-ray component. Interestingly, a component 
with large U, consistent with what we find in our model, is also 
invoked in th e NLR to produce optical high -ionization lines 
jKorista & Ferland.1989: Kraemer et al-lOm . 

An easy solution would be that all the components with 
different ionization parameters are spatially separated, possibly 
radially distributed, but the fact that the [O III] to soft X-ray ra- 
tio presented in this paper remains constant with the distance 
clearly favors the coexistence of various ionization parameters 
at each radius. This would imply a range in density and tem- 
perature at any given location, as already suggested, for exam- 



ple, bv lKinkhabwala et alJ lEooi andlO gle et all (l2003l) . Such 
a distribution m ay naturally arise as a r esult of thermal instabil- 
ities in outflows llKrolik&Krissl200 l'). In any case, the clumpy 
appearance of the [O III] emission regions is very clear from the 
HST images (see Fig.[2and|2}, thus implying that our model in 
terms of a single cone of gas with density smoothly decreasing 
outwards is oversimplified. Again, the role of radio jets in con- 
fining different components may be crucial. Unfortunately, the 
limited Chandra spatial resolution does not allow us to com- 
pare the substructures of the soft X-ray emission with those of 
the NLR and the radio jets, thus preventing us from drawing 
any firm conclusion on this issue. 

6. Conclusions 

We have pre sented a sample of e ight Seyfert 2 galaxies in- 
cluded in the lSchmitt et a l. ( 2003a) catalog of extended [O III] 
images, with a Chandra observation. All the sources but one 
present extended emission in soft X-rays, whose dimension and 
shape closely matches that of the [Olll]. The only exception, 
NGC 4507, does not show clear evidence for extended X-ray 
emission, but it is still consistent with the observed [O III] emis- 
sion. 

The spectral analysis of the sources suggests that the most 
likely origin for the soft X-ray emission is in a gas photoion- 
ized by the nuclear continuum. The clearest piece of evidence 
comes from the 190 ks combined RGS spectrum of Mrk 3, 
which unambiguously appears to be produced in a photoion- 
ized gas with an important contribution from resonant scatter- 
ing. In the other cases, where only spectra with CCD resolution 
are available, a 'scattering' model (a powerlaw plus emission 
lines) is to be preferred to a 'thermal' model either on statisti- 
cal grounds or because of unphysical best fit parameters of the 
latter (quasi-zero abundances). 

Therefore, we tested with the code CLOUDY if a solution 
in terms of a gas photoionized by the nuclear continuum which 
produces both the soft X-ray and the [O III] emission is tenable. 
Such solutions exist and require the density to decrease with 
radius. In particular, the fact that the observed [O III] to soft 
X-ray ratio is fairly constant up to large radii suggests that the 
density law should be close to r^^, similarly to what found in 
the NLR of some objects. 

In any case, these models are clearly oversimplified and do 
not take into account a number of details. In particular, there 
is evidence, both in the optical and in the X-rays, that var- 
ious components are present at any radius, as also indicated 
by the dumpiness of the extended region when seen with the 
high spatial resolution of HST. In this sense, the role of ra- 
dio jets, almost always found to be very well correlated with 
the NLR/soft X-ray emission, is still under debate and could be 
very important at least for the morphology of the gas. However, 
the investigation of more complex models would require better 
data quality, which will unfortunately not be available for many 
years to come. 
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Fig. 3. Chandra data and best fit models for the sources in our sample analysed for the first time in this paper. The imier box 
shows the 'local' fit for the iron Une, along with a dotted line to show the energy of the neutral iron Ka Une at the redshift of the 
source. 
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Fig. 4. Mrk 3: XMM-Newton 190 ks combined RGS spectrum, plotted in the rest frame of the source (12-25 A). The strongest 
emission features are labelled. 



